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A model consisting of a mixture of superconducting and quantum links is proposed to describe the 
integer quantum Hall transition. The quantum links correspond to tunneling of electrons between 
trajectories trapped in adjacent potential valleys, while the superconducting links correspond to 
merging of these trajectories once the Fermi energy crosses the saddle point energy separating 
the two valleys. The quantum Hall transition in this model corresponds to percolation of the 
superconducting links. Numerical calculations and scaling analysis using two different approaches 
yield the critical exponent v w 2.4 and a two-peak conductance distribution at the critical point. The 
role of quantum coherence is discussed, and an explanation of experimental observations claiming 
different universality class for the quantum Hall transition is suggested. The model suggests that 
the critical behavior of the superconductor-insulator transition from the insulating side is in the 
same universality class as the quantum Hall transition. 

PACS numbers: 71.30.+h,73.43.-f,73.43.Nq,74.20.Mn 
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The integer quantum Hall (QH) effect is one of the 
most studied examples of a second order quantum phase 
transition [l|. As such it is characterized by a diverging 
length scale £, describing the typical decay of a wavefunc- 
tion, or of the conductance. The physics describing the 
transition seems well understood, and most numerical in- 
vestigations in various models yield v ~ 2.4 0], where v 
describes the divergence of £ at the critical point. Experi- 
ments, however, disagree on the critical behavior. While 
some indeed yield an exponent close to the theoretical 
prediction |3|, others yield an exponent v ~ 1.3, close to 
the classical percolation critical exponent, v p = 4/3 [4|. 
In fact, the role of percolation in the QH transition has 
been experimentally demonstrated [5j- Yet other experi- 
ments claim that there is no critical behavior |g|. 

Similar lack of clarity exists with respect to 
the critical exponent describing the two-dimensional 
superconducting-insulator transition (SIT), which can be 
induced either by a continuous change of him thickness, 
or by a magnetic field. In both cases some experiments 
report an exponent v close to 1.3 on the insulating side 
|7j, while others yield v ~ 2.8 |8j. On the superconduct- 
ing side only v ~ 1.4 has been reported. Again, similarly 
to the QH case, some experiments give evidence of no 
critical behavior jj|] • The relevance of the percolation of 
superconducting areas to the SIT has been pointed out 
long time ago JlOj, and has been demonstrated experi- 
mentally both in granular systems [llj and in amorphous 
systems |l2l. More recent scanning tunneling measure- 
ments [l3j have directly demonstrated the separation of 
these systems into superconducting areas and metallic or 
insulating ones. 

In this work we develop a generic model to describe the 
QH transition. This model, which is directly related to 



the SIT, and to quantum percolation, clarifies the inter- 
play of classical percolation and quantum tunneling and 
interference in describing the critical behavior of the QH 
transition. In addition to providing an estimate of the 
critical exponent and the fixed-point conductance distri- 
bution at the critical point, it enables us to investigate 
the role of tunneling and interference on the critical be- 
havior, suggesting an explanation for the appearance of 
the classical critical exponent in the QH transition in 
some experiments Q. Moreover, the model suggests that 
the SIT, on the insulating side, belongs to the same uni- 
versality class as the QH transition, explaining the close- 
ness of the numerical and experimental estimates of the 
critical exponents of both transitions. 

The model proposed here is based on the following pic- 
ture. In strong magnetic fields electrons perform small 
oscillations around equipotcntial lines. At low energies, 
their trajectories are trapped around potential valleys, 
with quantum tunneling occurring between adjacent val- 
leys. We associate each such potential valley with a dis- 
crete lattice site, where nearest neighbors are connected 
by links that correspond to random tunneling between 
these trapped trajectories (see Fig. 1). As the Fermi en- 
ergy e_F rises and crosses the saddle point energy separat- 
ing two neighboring valleys, the two isolated trajectories 
coalesce, the electron can freely move from one valley to 
its neighbor and the link connecting them becomes per- 
fect or "superconducting" (SC). The QH transition oc- 
curs when one trajectory spans the whole system, which 
in this lattice model corresponds to percolation of the SC 
links. The critical behavior, however, will be determined 
by the quantum mechanical processes involved in trans- 
port of an electron from one side of the system to the 
other, namely tunneling and interference. 



For simplicity, the valleys are mapped onto sites on a 
square lattice. The saddle-points, which are mapped onto 
links in the lattice, have random energies e,-, sampled 
from some distribution G(ej). For a given Fermi energy, 
the transmission through the link for e F < £j is given by 



T(e F ) = exp [-a(ei - e F )} 



(1) 



describing tunneling through a parabolic barrier, with a 
a positive constant. For e F > ej, on the other hand, the 
transmission through the link is perfect (i.e. T(e F ) = 1) 
and the link is considered to be SC. At low Fermi energy, 
very few SC links exist, and their concentration p(e F ), 
given by p(e F ) = J_ F G(e)de, increases with e F . When 
p(e F ) = p c = 4f) the critical threshold for the square 
lattice, the SC links percolate through the system and the 
localization length diverges, signaling the QH transition. 
We will be interested in the critical behavior near this 
phase transition point |14|. 




FIG. 1: Mapping of the quantum Hall system onto a dis- 
crete model, (a) Equipotential lines for five different Fermi 
energies. At low Fermi energies electron trajectories are con- 
fined to potential valleys. Each such valley is denoted by a 
point in the discrete model (b)-(f). In order to move from one 
such valley to another, the electron has to tunnel through a 
saddle point - dotted links in (b)-(f). As the Fermi energy 
rises above the saddle point energy these valleys are joined 
and electron can move without dissipation from one valley to 
another - solid links in (b)-(f). Figures (b)-(f) correspond 
to a system with Fermi energy up to the equipotential line 
pointed by the appropriate letter in (a). Note that in (d) the 
solid links percolate, meaning that at this Fermi energy an 
electron can traverse the whole system on an equipotential 
line, corresponding to the quantum Hall transition. 

Two computational approaches for solving this quan- 
tum mechanical problem are presented. The first one 
employs a scattering matrix formalism. Each link in the 
lattice carries two edge states from neighboring valleys, 
moving in opposite directions (see Fig. 2a). The scat- 
tering matrix for each link is characterized by complex 
transmission and reflection amplitudes associated with 



electron tunneling between adjacent valleys. The tunnel- 
ing probability is determined as follows: for each link, a 
saddle-point energy ti is randomly chosen from a uniform 
distribution G(ei) on [— ^-, ^-] (so that the critical en- 
ergy is e c = 0, corresponding to p c = ^-). If e F > e^ then 
T = 1, otherwise the transmission probability through 
the link is given by Q ■ The allowed phases of the matrix 
elements are chosen randomly from a uniform distribu- 
tion between and 2tt. 

For perfect links (T = 1), the electron follows these 
edge states from one valley to another, according to their 
chirality. It is easy to see (Fig. 2b) that when these per- 
colate, we have two edge states propagating through the 
system in opposite directions, without scattering between 
left going and right going channels. Consequently, the 
percolation point corresponds to the quantum Hall tran- 
sition. 



X. 


t 


" 


;■ 


■' 






■ 


' 


• 






■ 


■ . 


- . 


X 


- 


. 


. 


'- 



K 



V x] 



3D 



IX 



DDDQ 



(b) 



FIG. 2: The scattering matrix approach: each link carries 
two counter propagating edge modes (a). A non zero trans- 
mission (broken lines) allows electrons to tunnel between adja- 
cent sites (potential valleys). When the transmission is unity 
(bold lines in (a)), these two valleys merge, and an edge state 
can freely propagate from one to another. A percolation of 
these perfect transmission links (b) correspond to an edge 
state propagating through the system without backscatter- 
ing. 

We evaluate the transmission T(L) = exp((logT)), 
through a system of the geometry depicted in Fig. 2, of 
linear sizes L = 10, 15, ... , 50, averaged over 5000 real- 
izations, for different values of the Fermi energy near the 
critical point. Fig. 3 displays the raw data. In (a) we plot 
the dependence of the transmission on Fermi energy (in 
terms of probability of SC links) for different sizes. The 
existence of a critical point, where the curves cross and 
the transmission is independent of length, is evident. In 
(b) we show the dependence of the transmission on length 
for different concentrations. The inset depicts the log of 
the transmission as a function of length. The straight 
lines demonstrate that indeed, T(L) oc exp(— L/^(e F )). 
All these curves coalesce (Fig. 4(a)), after scaling the 
system length the energy-dependent localization length 
£( e X- £( e X is found to diverge at the critical point with 
an exponent 2.34 ± 0.1, in agreement with previous nu- 
merical approaches. The inset of Fig. 4(a) displays the 
distribution of transmissions at the critical point. The 
double peak structure agrees w ith p revious numerical cal- 
culations for the integer QHE 
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FIG. 3: Numerical results of the scattering matrix approach 
for a = 12 (see Eq. 0. (a) The transmission coefficient T as 
a function of Fermi energy for system lengths L = 10, . . . , 50. 
Inset: The dependence of the critical point on the parameter 
a . Bold points are raw data and the solid line is the analytic 
curve (see text), (b) The transmission coefficient T as a func- 
tion of L for different concentrations, p — 0.41,0.4, . . . ,0.3. 
Inset: the same data on a semi-log plot demonstrating an 
exp(— L/£(eF)) dependence of the transmission. 
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FIG. 4: (a) The raw data curves for T (Fig. 3) are collapsed 
onto a single curve by scaling each curve by an energy depen- 
dent localization length. This localization length diverges at 
the critical energy with an exponent 2.34±0.10. The fixed dis- 
tribution of T at the critical point is depicted in the inset, (b) 
Scaling of the data obtained by the tight-binding approach. 
For energies e = —0.050, —0.055, . . . , —0.1 the data collapse 
yields the exponent v = 2.43 ± 0.1. Inset: the tight-binding 
approach - sites connected by perfect links are merged and 
then the transmission coefficient of the modified lattice is cal- 
culated using the tight-binding hamiltonian @. (c) Data 
collapse with the same exponent as in (b), including a larger 
range of energies e = —0.01, —0.02, . . . , —0.2. This data can 
be better fitted with an exponent v = 4/3 (d). 

The exponent v is found to be independent of the value 
of a in l|I]l. The critical probability, on the other hand, 
does depend on a (inset of Fig. 3(a)). This dependence 
can be understood as follows: The transmission through 
the system at the critical point, Tq, is independent of 
length and of the value of a and is given by an aver- 
age over the fixed distribution (inset of Fig. 4(a)). As 
the transmission rises with Fermi energy ej?, the critical 
Fermi energy for a given a, e c (a), is determined by the 
value of 6f at which the transmission becomes equal to 
Tq. Since the distribution of link transmissions is expo- 



nential, the conductance of the whole network is deter- 
mined by the critical transmission, the smallest trans- 
mission value such that the links with better transmis- 
sion percolate through the system [lfj. This transmis- 
sion is given by exp[a(e c (a) — e c )], with e c = 0, the 
classical critical Fermi energy. Equating this transmis- 
sion to To, one immediately finds e c (a) — log(Tb)/a, or 
p c {a) = 4r + log(To)/a. This function is plotted in the 
inset of Fig. 3(a) (solid curve), with T = 0.4, and demon- 
strates an excellent fit to the data. The average value 
of log(T) over the fixed distribution found above gives 
To = 0.32, in rough agreement with the argument. 

In order to study a model where the critical point is in- 
dependent of parameters, we develop a second approach, 
employing a tight-binding Hamiltonian 



n 



E 
<hj> 



■Jijcjcj + H.c. 



(2) 



Determination of the parameters vy follows the pro- 
cedure described above for determining the scattering 
matrix parameters. We associate a random saddle-point 
energy e^ with each link joining the lattice sites i and 
j. When tij > ej? the hopping matrix element Vij for 
that link is chosen such that it will give the transmission 
coefficient Q . When e^ < ep the link is considered per- 
fect, or SC. Since no value of the hopping matrix clement 
corresponds to perfect transmission in the tight-binding 
approach, then when two sites are joined by such a SC 
link, they are merged into a single site. This guaran- 
tees that only when the SC links percolate, the conduc- 
tance of the system will not decay with length, so that 
p c will be independent of a. The conductance is now cal- 
culated numerically, by attaching a single-channel wire 
to both sides of the lattice, and calculating the trans- 
mission through the system for an electron at the Fermi 
energy. This calculation involves two stages (see inset of 
Fig. 4(b)): (a) identifying clusters of sites connected by 
SC links as a single site, and (b) calculating T by solv- 
ing the tight-binding equations in the reduced space [TtJ • 
Data is generated for lattice length L = 6, 10, ... , 34, 
and averaged over 1500 realizations. Following the same 
scaling procedure as described above, the curves collapse 
onto a single curve, a procedure that yields the critical 
exponent v = 2.43 ±0.1 (Fig. 4(b)). We also calculate 
the conductance distribution at the critical point and find 
that it shows a two-peak structure, similar to that of the 
S-matrix approach. Both the critical exponent and the 
critical point, p c = ^- are independent of the parameter 
a, for a large range of values [18|. We also check the 
dependence of the critical exponent v for different func- 
tional forms of the tunneling probability T(e) in addition 
to that described by 0). Using, for example , th e func- 
tional form obtained by Fertig and Halpcrin ljj], yields 
very similar results, with a value v = 2.49±0.1. A similar 
value is obtained when the link transmission coefficients 



are taken from a uniform distribution. Thus, the model 
does not depend on the exact functional form of T(e), or 
on its distribution. The change in the critical exponent 
from its classical counterpart can be traced to the fact 
that as the Fermi energy, or concentration of SC links, 
change, the transmission amplitude through the quantum 
insulator is modified, as the bottle-neck link changes. We 
also note that repeating the same calculations without 
random phases does not yield any reasonable scaling of 
the data. Thus interference effects are crucial in obtain- 
ing the correct critical behavior. 

We also note that if one tries to employ the scaling 
analysis described above including energies that are not 
too close to the critical one, the collapse of the data be- 
comes worse and in fact can be better fitted with the 
classical exponent v — 4/3 (Fig. 4(c) and (d)). This is 
consistent with the above picture as the difference be- 
tween the classical and the quantum exponents arise due 
to the additional dependence of the tunneling amplitude 
on e — e c , a dependence that becomes less relevant away 
from the critical point. This observation may also explain 
why some experiments report a classical percolation ex- 
ponent for the QH transition [4J. 

Since this model is practically identical to a model de- 
scribing the SIT, it is not surprising that the observed 
critical exponents (on the insulating side) for the latter 
are very similar to those observed for the former. How- 
ever, unlike the QH transition, the SIT is not symmetric 
around the critical point, which is closer to the model de- 
scribed here. In fact, in the superconducting phase the 
mechanism described above for the change in the critical 
exponent from its classical percolation value due to the 
change in the transport in the insulating phase does not 
have any effect. This is consistent with the experimen- 
tal observations of a critical exponent close to v w 1.4, 
even for samples where a larger exponent was observed on 
the insulating side [SJ- In addition, this relation between 
these two apparently different transitions also explains 
the fact that the critical resistance observed at the SIT 
seems to be distributed around values somewhat larger 
than /i/4e 2 , expected from duality[2fjj. This relation can 
be further investigated experimentally by studying the 
full resistance distribution at the SI threshold, similarly 
to such experiments for the QH transition [2 1|. The simi- 
larity of the QH transition and the SIT has already been 
pointed out in |22j, who emphasized the non-critical be- 
havior and the role of decoherence. In fact, the model 
introduced in this work allows the introduction of de- 
phasing in a straightforward way, by attaching current- 
conserving, phase-breaking reservoirs to some fraction 
of the links |23j . The interplay between dephasing and 
classical-quantum crossover will be explored in a future 
communication J24| . 
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